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The Jahn-Teller distortion, by its very nature, is often at the heart of the various electronic 
properties displayed by perovskites and related materials. Despite the Jahn-Teller mode being non¬ 
polar in nature, we devise and demonstrate in the present letter an electric field control of Jahn-Teller 
distortions in bulk perovskites. The electric held control is enabled through an anharmonic lattice 
mode coupling between the Jahn-Teller distortion and a polar mode. We conhrm this coupling, 
and explicitly an electric held effect, through hrst principles calculations. The coupling will always 
exist within the Pb2im space group, which is found to be the favoured ground state for various 
perovskites under sufficient tensile epitaxial strain. Intriguingly, the calculations reveal that this 
mechanism is not only restricted to Jahn-Teller active systems, promising a general route to tune 
or induce novel electronic functionality in perovskites as a whole. 


Perovskites, and related materials, are fascinating sys¬ 
tems exhibiting a diverse collection of properties, includ¬ 
ing ferroelectricity, magnetism, orbital-ordering, metal- 
insulator phase transitions, superconductivity and ther- 
moelectricityi. Despite the wide range of physical be¬ 
haviour, a common point at the origin of many of them 
can be identified as being the Jahn-Teller distortion^i^. 
The Jahn-Teller distortion is itself intimately linked to 
electronic degrees of freedom, since traditionally it mani¬ 
fests to remove an electronic degeneracy, opening a band 
gap and favouring a particular orbital ordering, which in 
turn can affect magnetic ordering. Furthermore it plays 
an important role, for example, in colossal magnetoresis¬ 
tance phenomena in doped manganites — , superconduc¬ 
tivity^!^ or the strong electronic correlation observed in 
the thermoelectric NaCo 02 family^. 

It would be highly desirable, for device functionality 
for example, to be able to tune the Jahn-Teller distor¬ 
tion and hence its corresponding electronic properties, 
with the application of an external electric held. How¬ 
ever, Jahn-Teller distortions are non-polar and hence not 
directly tunable with an electric held. One strategy to 
achieve an electric held control of Jahn-Teller distortions 
might begin by engineering a polar structure. Usually, 
perovskites adopt a non-polar Pbnm ground state, re¬ 
sulting in a combination of three rotations of the oxygen 
octahedra {a~a~c^ in Glazer’s notation^), also called an- 
tiferrodistortive (AFD) motions. Unfortunately, these 
AFD distortions are known to often prevenli^ the ap¬ 
pearance of the polarization in the material. The use 
of anharmonic lattice mode couplings between polar and 
non-polar lattice distortions is a promising pathway to 
engineer “improper” ferroelectricity in layered perovskite 
derivativesiSrJ^, and furthermore to achieve novel func¬ 
tional possibilities via these couplingsiir— . Interest¬ 
ing anharmonic couplings do not only exist in layered 
perovskites, but can also appear in bulk ABO 3 per- 
ovskite a^^i^^i^^i^^~ — 

Following this spirit, achieving an electric held control 
of Jahn-Teller distortions necessarily requires the identi- 
hcation of a material exhibiting a relevant coupling be¬ 



FIG. 1. Schematic view of the four lattice distortions in¬ 
volved in the Pb2im phase of perovskites under tensile epitax¬ 
ial strain, a) Polar distortion (irreps T^) b) Qj Jahn-Teller 
distortion (irreps ) c) Anti-polar A distortion (irreps ) 
d) a^a^c^ (jif antiferrodistortive motion (irreps Mj). 


tween the polarization and the Jahn-Teller distortion, 
which to the best of our knowledge has not yet been 
discovered in bulk perovskites^^. In the present letter we 
identify such conditions, and demonstrate explicitly an 
electric held control, in bulk perovskites using a combi¬ 
nation of symmetry analysis and hrst principles calcula¬ 
tions. 


RESULTS 

The two required lattice distortions are pictured in hg- 
ure [TJ a) and b). From a symmetry analysis of the dif¬ 
ferent tilt patterns in combination with possible Jahn- 
Teller lattice distortions, Howard and Carpenter pointed 
out that a Jahn-Teller distortion pattern automatically 
appears when considering the Pbnm symmetry^^, later 
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explained in terms of anharmonic couplings with AFD 
motionsi^. As a consequence, a Jahn-Teller distortion is 
not necessarily electronically driven, but can instead arise 
from lattice mode couplings in which case a splitting of 
the electronic states may develop even in the absence of 
a degeneracy of states. The present Jahn-Teller lattice 
motion corresponds to a Q 2 mode as defined by Goode- 
nougb^, corresponding to two B-0 bond length contrac¬ 
tions and two B-0 elongations. This motion orders at 
the M point of the Brillouin zone and hence consecutive 
layers along the c axis of the Pbnm phase present in- 
phase distortions. Consequently, this motion is labelled 
QJ (irreps throughout the whole manuscript. 

Starting from the ideal Pm3m cubic perovskite phase, 
the condensation of the polar mode P (irreps T^) and the 
JT mode lowers the symmetry to a Pb2im phase, a 
polar subgroup of Pbnm. We then perform a free energy 
expansion^ in terms of the possible lattice distortions in 
this new phase and we identify, among all the possible 
terms, some intriguing couplings: 

J^(xPQ+A + p2Q+<^++p^+yl + Q+0+yl2 (1) 

In this phase, the first two terms of equation [T] provide a 
link between the polarization and the Jahn-Teller distor¬ 
tion. These terms also involve two additional distortions: 
one anti-polar A motion pictured in figure [1] c) and one 
a°a°c+ AFD motion (labelled 0+) pictured in figure [T] 
d). Among all the terms, the lowest order trilinear term 
of the form provides the desired direct coupling 

between the polarization and the JT distortion. Thus, 
acting on the polarization with an external electric field 
may modify the amplitude of the JT motion, and there¬ 
fore all related electronic properties. 

However, as previously discussed, the Pb2im symme¬ 
try is not the common ground state in bulk perovskites^l. 
Strain engineering, through thin film epitaxy for exam¬ 
ple, can provide a powerful tool to unlock a polar mode in 
perovskite o^^i^°~ — . This is the case for BiFeOs that was 
recently proposed to adopt an unusual Pb2im symmetry 
under large epitaxial tensile strai n^^i^^i^^ . This particular 
phase was shown to develop both polar, anti-polar and 
a°a°c“'' AFD motions^^, which were later demonstrated 
to be coupled together through the third term of eq. 
Amazingly, the authors reported the existence of an or¬ 
bital ordering of the Fe^"*" 3d orbitals, explained from 
the coexistence of the polar and the anti-polar motion 
yielding a particular lattice distortion pattern^S. This 
orbital-ordering is unusual since in this system no Jahn- 
Teller effect is required to form a Mott insulating state 
(Fe^’*' are in a half hlled - high spin conhguration). 

A Jahn-Teller effect or distortion are yet to be reported in 
the Pb2im phase of BiFeOs to the best of our knowledge. 
From our symmetry analysis, we clearly demonstrate 
that as this Pb2im develops the three aforementioned 
distortions (P, A and the free energy of eq. [T]is au¬ 
tomatically lowered through the appearance of a fourth 
lattice distortion: a Jahn-Teller Qj motion. Therefore, 
whilst it may not itself be unstable, the Jahn-Teller mo¬ 


tion is forced into the system via this “improper” mech¬ 
anism arising from the trilinear couplingiS. This result 
clarifies the origin of the unusual orbital-ordering dis¬ 
played by BiFeOa and moreover, it provides a pathway 
to achieve an electric field control of the orbital-ordering 
in bulk perovskites. 

The predicted highly strained Pb2im phase in bulk 
perovskites is not restricted to BiFeOa, and it was pre¬ 
dicted to occur also in some titanates (CaTiOa and 
EuTiOa)^, in BaMnOg^ and even in a Jahn-Teller ac¬ 
tive compound TbMnOg^. The highly strained bulk per¬ 
ovskites are then an ideal playground to demonstrate our 
coupling between the polarization and the Jahn-Teller 
distortion. In order to check the generality of our con¬ 
cept, we propose in this letter to investigate several types 
of highly strained perovskites on the basis of first princi¬ 
ples calculations: i) non magnetic (NM) SrTiOa (t^gCg); 
ii) magnetic BaMnOg^ ^-ud BiFeOs iii) 

Jahn-Teller active YMnOs (t^gCg)- 

We begin by investigating the possibility of a Pb2im 
ground state under large epitaxial tensile strain (the 
growth direction is along the [001] axis of the Pbnm 
structure). Beyond around 5% tensile strain, the four 
compounds indeed develop the desired Pb2im ground 
state. Strained BaMnOa (ferromagnetic FM) and 
YMnOa (G-type antiferromagnetic AFMG) exhibit a 
different magnetic ground state compared to the bulk 
(AFMG and E-type antiferromagnetic - ttii zig-zag 
chains coupled antiferromagnetically along the c axis 
- respectively) while BiFeOa (G-type antiferromagnetic 
AFMG) remains in its bulk magnetic ground state. We 
then perform a symmetry mode analysis with respect to a 
hypothetical PA/mmm phase (corresponding to Pm3m 
for unstrained bulk compounds) in order to extract the 
amplitude of the relevant lattice distortions^ (see ta¬ 
ble |T|. As expected, the four materials develop the re¬ 
quired distortions, and amazingly, the magnitude of the 
QJ Jahn-Teller distortion is relatively large, being for 
instance of the same order of magnitude as the one de¬ 
veloped in the prototypical Jahn-Teller system LaMnOa 
(around 0.265 A^i). Additionally, the values of the 
spontaneous polarization are rather large, reaching 76 


SrTiOs BaMnOs BiFeOs YMnOs 


strain 

(%) 

-h7.35^ 


-15.8^ 

-b4.0^ 

magnetism 


NM 

FM 

AFMG 

AFMG 

P (Ps") 

(A) 

0.615 

0.421 

0.346 

0.753 

{jiC.cm ^) 

76 

45 

29 

7- 

(M^) 

(A) 

0.232 

0.190 

0.644 

0.737 

A(M+) 

(A) 

0.558 

0.217 

1.072 

0.940 

(M+) 

(A) 

0.640 

0.059 

1.668 

1.733 

gap 

(eV) 

3.02 

0.28 

1.88 

1.88 


TABLE I. Epitaxial strain (%), magnetic ground state, am¬ 
plitudes of distortions (A) and electronic band gap value (eV) 
for each material. We emphasize that only the relevant dis¬ 
tortions are summarized in the present table. 
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FIG. 2. Potentials with respect to the amplitude of dis¬ 
tortions of the four lattice motions producing the required 
Pb2im for SrTiOs (black filled circles), BaMnOs (red filled 
squares), YMnOs (blue filled diamonds) and BiFeOs (green 
filled triangles) starting from the ideal P4/mmm phase. 



FIG. 3. Projected density of states on the d levels on two 
neighboring B sites in the (a6)-plane of SrTiOa, BaMnOa, 
BiFeOa and YMnOa. Local axes of the orbitals are displayed 
on figure [T] The Fermi level is located at 0 eV. 


for SrTiOa for instance. Despite being highly 
strained, all materials remain insulating, adopting rea¬ 
sonable electronic band gap values (see table H]). 

To shed more lights on the origin of this unusual 
Pb2i'm phase we compute the energy potentials with 
respect to the four distortions by condensing individu¬ 
ally each modes in an hypothetical PA/mmm phase (see 
figure [Hi. Surprisingly, the appearance of the Pb2im 
phase is rather different for the four materials. SrTiOs 
and BaMnOs only exhibit a polar instability, producing 
an Amm2 symmetry, consistent with previous reports 
of a polar phase for these two materials under tensile 
strain^Sii^. Computing the phonons in this particular 
Amm2 symmetry, only one hybrid unstable phonon mode 
is identified for these two materials, having a mixed char¬ 


acter between the A, 0+ and Qj distortions. For BiFeOs 
and YMnOs, the a^a^c'^ AFD motion is already unsta¬ 
ble, which is expected since the Pb2im symmetry for 
these two systems is derived from their bulk R3c/Pbnm 
phases^. Additionally, the JT lattice distortion is also 
unstable in the PA/mmm phase of YMnOs and ap¬ 
pears as an electronic instability^, which is expected 
since YMnOs is known to be Jahn-Teller active in the 
bulk. We emphasize at this stage that the polar mode 
in BiFeOs (and YMnOs) is not unstable and therefore 
highly strained BiFeOs appears as an improper ferro¬ 
electric in contradiction to reference [13. Computing the 
phonons in the intermediate strained Pbnm phase of both 
BiFeOs and YMnOs compounds reveals only one hybrid 
unstable mode, having a mixed character between P and 
A distortions.Despite the apparent universal stability of 
this highly strained polar phase, the mechanism yielding 
it is suprisingly different between the compounds. 

Regarding the electronic structure, we checked for 
the appearance of an orbital ordering as observed in 
BiFeOs^^. For the four compounds we report the pro¬ 
jected density of states on the d levels of two neighboring 
B sites in the (a6)-plane (see figure |3]). For SrTiOs, a 
splitting of the t 2 g states, and especially between the dxz 
and dyz orbitals, located at the bottom of the conduction 
band arises. For BaMnOs and BiFeOs, a similar splitting 
between the t 2 g levels is observed near the Fermi level, 
even if it is less pronounced for BaMnOs since it has the 
smallest Q 2 distortion. Finally, YMnOs displays an or¬ 
bital ordering of the Cg levels with predominantly d^^-y^ 
occupation. This splitting is known to result of the Jahn- 
Teller distortion in this A^+Mn^+Os class of material^. 
Additionally, an orbital ordering of the Og levels is occur¬ 
ring both in the conduction and the valence bands. To 
prove that the Jahn-Teller distortion, and not another 
motion, is solely responsible for the orbital ordering we 
have condensed all the modes individually and studied 
the density of states (see supplementary figure 1). 

Finally, we explicitly demonstrate the electric field con¬ 
trol of Jahn-Teller distortion by computing, as a function 
of the electric field E applied along the direction of the 
spontaneous polarization, the evolution of both the JT 
distortion and the band gap. Results are displayed in 
figure m The Jahn-Teller distortion is effectively altered 
by the application of an electric field along the polar axis 
through the first and second terms of equation [T] As the 
electric field increases, the amplitude of the JT distor¬ 
tion is either amplified or decreased, being renormalized 
to around 175% for SrTiOs for an electric field around 
20 MV.cm~^. The largest effect is however reached for 
YMnOs which displays a renormalization of 130% un¬ 
der moderate electric field (around 5 MV.cin”^). There¬ 
fore, this renormalization of the JT distortion has con¬ 
sequences for instance on the electronic band gap value, 
with an opening/closure around 0.6 eV for YMnOs or 
0.25 eV for SrTiOs. It is then possible, through the 
coupling between the polarization and the Jahn-Teller 
distortion to act on the electronic band gap, and one 
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FIG. 4. Relative electric field effect on the amplitude of the 
Jahn-Teller distortion (top panel) and the electronic gap value 
(bottom panel) on the four different compound. 

may imagine a possible control of a Metal-Insulator phase 
transition, optical properties, or photovoltaic efficiency, 
for instance. 

In conclusion, we have demonstrated in the highly 
strained Pb2im phase of bulk perovskites the existence 
of a coupling between a polar mode and the Jahn-Teller 
distortion. This improper anharmonic coupling, estab¬ 
lished on universal symmetry arguments, enables an elec¬ 
tric field control of the Jahn-Teller distortion, even in the 
case of non electronically Jahn-Teller active systems. The 
generic mechanism may open novel functionalities in per¬ 
ovskites as it will have consequences on related electronic 
properties as proposed in the present letter. For instance, 
such couplings may allow the tuning of Metal-Insulator 
phase transitions, the control of orbital orderings, optical 
properties, or of the photovoltaic efficiency. 


METHODS 

First-principles calculations were performed with the 
VASP packag o^^ds . yvb used the PBEsol^-l-U frame¬ 
work as implemented by Lichtenstein’s method^ (see the 
supplementary material for a discussion on the choice 
of the U and J parameters). The plane wave cut-off 
was set to 500 eV and we used a 6x6x4 k-point mesh 
for the 20 atom Pb2im phase. PAW pseudopotentials^ 
were used in the calculations with the following valence 
electron configuration: 3s^3p®4s^ (Sr), 4s^4p®5s^ (Ba), 
4s^4p®5s^4d^ (Y), 6s^6p^ (Bi), 3p®4s^3d^ (Ti), 3p®4s^3(i® 
(Mn), 3p®4s^3(i® (Fe) and 2s^2p^ (O). Spontaneous po¬ 
larizations were computed using the Berry-phase ap¬ 
proach and phonons and Born effective charges were com¬ 
puted using the density functional perturbation theory^. 
The electric field effect was modelled using a linear re¬ 
sponse approach by freezing-in some lattice distortion 
into the syste m^^d^ . Symmetry mode analyses were per¬ 
formed using the Amplimodes software from the Bilbao 
Crystallographic serve r^^d^ . 
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